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Abstract— Electrocardiogram (ECG) is a vital signal which 
represents the state of health of astronauts in a manned 
spaceflight mission and hence must be continuously acquired and 
transmitted throughout the mission. Telemetry bandwidth is a 
premium resource in such applications. Traditional method of 
sampling at Nyquist rate is highly bandwidth inefficient. Hence it 
is advantageous to use Compressive sensing (CS) technique to 
optimize data at measurement point itself. Conventional CS 
techniques employ computationally intensive measurement 
matrices which are not hardware efficient for both acquisition as 
well as recovery. In this work, a hardware efficient scheme with 
use of a sparse binary measurement matrix is proposed. The 
signal is recovered from the compressive measurement using a 
constraint function based on inverted Laplace distribution 
function. Gradient decent method is used to recursively recover 
the original ECG signal from the compressive measurements in 
an efficient manner. Apart from this, a projection scheme is also 
proposed to minimize the recovery error. The proposed scheme 
was extensively evaluated with different ECG samples with 
different compression ratios. Finally, the proposed scheme was 
benchmarked with Approximated LO norm based method and it 
is found to perform more efficiently in compressively sensing and 
recovery of ECG signals. 


Keywords— Electrocardiogram (ECG), compressed sensing 
(CS), discrete cosine transform (DCT), constraint function. 


I. INTRODUCTION 


The rigors of space travel during the ascent phase of the 
launch vehicle, in orbit phase weightlessness and the toughest 
re-entry phase back to earth puts considerable physiological 
stress on the astronauts in a manned spaceflight. Hence ECG 
which is a vital health parameter of astronauts are to be 
continuously monitored onboard and transmitted through 
telemetry to ground stations right throughout the mission. But 
the constraints of data rate necessitate use of compression 
techniques. The acquisition of signals at Nyquist rates and 
then applying conventional compression techniques is 
inefficient. Compressed sensing (CS) [1],[2] overcomes this 
limitation by directly acquiring data in a _ condensed 
representation. CS works with an important property of the 
input signal, called sparsity. If the samples generated by 
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sampling a signal is represented as a vector, and if the vector is 
having only K non-zero entries, then signal is called K-Sparse. 
CS theory doesn’t stipulate that signals be sparse in time- 
domain. Instead if a signal is sparse in a particular domain or 
basis, then CS technique can be applied in that basis to 
significantly reduce the sampling rate of the signal. 


Fei-Yun Wu et.al [3], proposed a compressive sensing 
technique for underwater acoustic signals, wherein a Gaussian 
random matrix optimized into a Gram matrix is used. The 
authors have proposed the use of Approximated Lo norm 
(ALO) for reconstruction of the signal. He Caiwei et.al [4] 
proposed an Information-Weighted Gaussian matrix (IWGM) 
to improve the quality of the recovered ECG signal from 
compressed sensing. Fahimeh Ansari-Ram and _ Saied 
Hosseini-Khayat [5] proposed ECG signal compression using 
compressed sensing with non-uniform binary matrices. The 
method employs transmitting the sample numbers of the QRS 
complex to the receiver for improved recovery of the region of 
interest. Zhimin Zhang et.al [6] proposed an ECG-CS scheme 
by mapping the sub-sampled signals to a 2D cut and align 
map. Matching Pursuit and block sparse Bayesian learning 
algorithms were used by them for recovery. These CS 
techniques published in literature cannot be directly used for 
the intended application of real time CS based ECG signal 
acquisition onboard space missions and its recovery in ground 
stations. The primary reason for this is that the published CS 
works [7]-[10] employ either a computationally intensive 
measurement matrices or use complex optimization techniques 
for these measurement matrices. Furthermore, the recovery of 
original signal from the compressively sensed signals requires 
iteratively longer and complex recovery algorithms. All of 
these put higher loads on the onboard hardware. This paper 
details an efficient CS based acquisition and reconstruction 
scheme for ECG signals. The scheme is presented in Section 
II. The experimental studies are presented in Section-III. 
Section-IV discusses the results and performance of the 
proposed scheme. 
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Il. THE PROPOSED SCHEME 


A. The Overall CS Framework 


Consider a discrete signal x, which is real-valued, finite 
length and is one dimensional in nature. That is, x €RN, and 
can be viewed as Nx1 column vector in 8" with elements x[n], 
where n= 1, 2, 3,...., N. Now suppose that a measurement 
system acquires the signal x having finite length N using M 
linear measurements, such that M<<N. This can be 
mathematically represented as: 


(1) 


where y has length M, ® is an MxN measurement matrix with 
sub-sampling rate R=M/N. It would certainly be felt that 
recovering the signal x with very less measurements is 
impossible given that the number of unknowns far exceeds the 
number of measurements. But, CS theory tells that it is indeed 
possible provided the signal x is sparse in some domain or 
basis. Such sparse representation may be obtained by some 
transform as represented by: 


Yux1 = Omxn Xnx1 


x=WYs (2) 
where is the NxN basis matrix of x and s is a Nx1 column 
vector of weighting coefficients such that 


c= (3) 


Here is the transpose of Y. From (3) it is clear that both x 
and s are representations of the same signal in two domains, x 
in say time domain and s in ¥ domain. Thus x is K-sparse in ‘¥ 
domain. The compressed sensing equation (3) can be thus re- 
written as: 


y = Ox = OVs = Os (4) 


where © = OY is anM x N matrix. 


The measurement process is assumed to be non-adaptive [9], 
which means ® is independent of the signal x. The number of 
required measurements M is given by [1]: 


(5) 


where C is a constant called observation factor usually set to 1 
and K is the number of non-zero elements in the signal 
representation. 

The central problem of CS as represented by (4) is thus a 
randomized dimensionality reduction problem. The challenges 
in CS can be categorized into two: 


N 
M = CK log (<) 


1. Design of a stable measurement matrix ® in such a 
way that the effect of the dimensionality reduction 
from xeRN to yeR™ on the information retrieval is 
nullified and 

2. Development of a reconstruction algorithm that can 


effectively recover x with N samples from just M 
measurement. 
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B. Proposed Measurement Matrix 


Random Gaussian matrices are widely used in CS 
applications as they satisfy Restricted Isometry Property (RIP) 
[1] with overwhelming probability and have minimal mutual 
coherence measure [1], which are necessary conditions for 
sparse signal recovery. However, limitations such as floating 
point representation, implementation complexity and storage 
make them practically unviable. Considering these aspects, 
sparse random binary matrices which are easier to represent 
and implement with comparable recovery properties with 
respect to random Gaussian matrices are proposed in this 
work. 

The sparsity degree of random binary measurement 
matrix, D=[bijJeRMN - P(dij=l)=p and P(ij=0)=/-p is 
controlled by the parameter p which defines the probability of 
an entry in the binary matrix to be one. For sparse binary 
matrices, the value of p is very small. An average of Np 
number of ones will be there in each row of the matrix. Thus, 
the proposed binary matrix based measurement system 
requires an average of M (Np-1) adders and no multipliers. It 
may be noted that a Gaussian matrix based measurement 
system requires M (N-1) adders and MN multipliers. It may 
also be noted that no complex optimization techniques usually 
employed for random Gaussian measurement matrices, is 
required in this case. Hence, this method is highly hardware 
friendly. Detailed simulations of the proposed method in 
MATLAB, using various ECG databases, show that this 
scheme works quite efficiently. 


C. Estimation and Recovery 


In order to achieve a refined recovery result of original 
ECG signal from its compressively sensed measurements, a 
suitable constraint function that closely approximates the LO 
norm is proposed in this work. The close approximation to LO 
norm is used since the LO norm based approach is Non- 
deterministic polynomial hard (NP hard). 

Normalized and inverted probability density function 
(pdf) of Laplace distribution of a random variable is proposed 
as a constraint function which will be used for the 
reconstruction of the sparse signals from the measurements. 

A random variable € has a Laplace distribution if its 
probability density function is: 


f= (Se, ite <u 


: (6) 
1\ -(+) 
f@=(S)e), te >a (7) 
where yw is a location parameter and 5 is a scaling parameter 
(with b>0). Using (6) and (7), the constraint function 
(Normalized and inverted Laplacian pdf) proposed in this work 

is given by: 


f(@®) =1 ele), ife<u (8) 


(9) 


e— 
b 
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The proposed constraint function defined by (8) and (9) for 
different scaling parameter b and location parameter uw = 0 is 
shown in Fig.1. 
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Fig. 1. Proposed constraint function 


Considering unavoidable noise in the overall system, the 
compressed sensing equation, y= Os as given in (4) is 
constrained by the equation: 


lly -Osll < B (10) 


where B represents the unavoidable noises. 


The reconstruction is based on searching for an optimal 
solution, 8 for the sparse representation s through convex 
optimization relation: 


f(@) such that |ly-—@s||< B (11) 


The solution for (11) is done in two parts. First, Gradient 
decent method [3] is used for finding the nearest local 
minimum for the function f (0) by computing the gradient of 
the function. As per the least square optimization technique, 
the recovered signal @ is given by: 


9=o0'y (12) 


where ©" is the Moore-Penrose pseudo-inverse defined by: 
o' = 0" (00")* (13) 


The k recursive iteration of the Gradient decent algorithm is 
defined by the equation: 


_ 0 
Ox+1 = 9% — O% 0 rq (F@)) (14) 


where A is the step size and o represents dot product. 
From (8), 


(FO) = fos) as) 
Substituting (15) in (14), we get: 


~sign(e)(°# 
O41 = Ox — orel sien(("5")) (16) 


In order to keep the estimation error under check, a 
projection method as in the work published by Fei-Yun Wu 
et.al [3] is also performed in parallel to this Gradient decent 
method. The projection is given by the relation: 


Oupdated = Ox41 + ty — P0x41) (17) 


The recovered ECG signal is obtained by taking the inverse 
discrete cosine transform (idct) by the equation: 


Z= idct(®updatea) (18) 


D. Overall Advantages of Proposed Scheme 


In comparison with existing CS methods, the proposed 
method gives improved overall recovery characteristics due to 
use of a constraint function that is a very close approximation 
to LO norm. Also, the use of sparse binary matrix yields a 
significant hardware advantage with regard to handling and 
storage. Due to these aspects, the proposed scheme is suitable 
for real time CS applications. 


Il]. EXPERIMENTAL STUDIES 


This section details the experimental studies carried out to 
evaluate the proposed CS framework for ECG signals. The 
proposed scheme was simulated in MATLAB R2014b 
environment with Intel i5-8250U CPU @1.6GHz and 8 GB 
RAM. 


ECG samples from MIT-BIH database was used as input 
signal to evaluate the performance of the proposed CS method. 
The ECG input sample consisted of 24 cycles of ECG 
waveforms with 10000 samples. The proposed method was 
coded in MATLAB and evaluated with the ECG database for 
different compression ratios. The compression ratio is defined 
as: 





CR= (- = ~) X 100 (19) 


where N is the finite length of the input signal and M is length 
of the compressed measurement. 


For simulation run for CR=50%, the typical settings for the 
variables in the simulation are N=400, M=200, A=0.01. The 
normalized recovery error is computed for every cycle of the 
ECG waveform. 








The robustness of the proposed scheme was also evaluated 
with anomalous ECG data like arrhythmia cases as well as 
ECG samples taken with motion artifacts. 


Authorized licensed use limited to: Auckland University of Technology. Downloaded on June 03,2020 at 15:46:00 UTC from IEEE Xplore. Restrictions apply. 


The DCT coefficients of the input ECG samples are shown 
in Fig.2. 
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Fig. 2. DCT of input ECG samples 


From Fig.2, it is quite evident that ECG signals are sparse 
in the DCT domain. The normalized recovery error for all the 
24 cycles was also computed to analyze the recovery error 
with different compression ratios. 


IV. RESULTS AND DISCUSSION 


This section details the results of the extensive simulations 
conducted to evaluate the performance of the proposed method. 
The proposed method was also benchmarked and compared 
with ALO method [3] and comparative results are brought out. 


The simulations were carried out with different ECG 
waveform samples including cases with arrhythmia as well as 
ECG acquired with motion artifacts. The simulation runs for 
each of these ECG samples were with different settings of 
compression ratios. 


The input ECG signal and recovered signal using the proposed 
scheme for a compression ratio of 50% is shown in Fig.3. 
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Fig. 3. ECG recovery with CR=50% 


The corresponding normalized recovery error is shown in 
Fig.4. 


0.012 T r 7 7 





0.01 4 


Magnitude 
o 
i=] 
8 
. 
1. 








0.002 











10 15 20 25 
Segment Number 


Fig. 4. Normalized recovery error. CR=50% 


The input ECG signal and recovered signal using the 
proposed scheme for a compression ratio of 60% and the 
corresponding normalized recovery error is shown in Fig.5 and 
Fig.6 respectively. 
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Fig. 5. ECG signal recovery with CR=60% 
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Fig. 6. Normalized recovery error. CR=60% 
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Table I : Comparative performance of proposed method 
with ALO method for ECG signal recovery 




















Compression Normalized Normalized 
ratio (CR) recovery error for | recovery error for 
ALO method [3] | proposed method 
25% 13x 10° 1.4x 104 
50% 0.04 8.5 x 107 
60% 0.4 0.09 
70% 0.85 0.44 
75% >1.2 0.65 

















ECG signal with arrhythmia condition and recovered signal 
using the proposed scheme for a compression ratio of 50% is 
shown in Fig7. 
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Fig. 7. ECG signal recovery arrhythmia case 
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Fig. 8. ECG signal recovery with motion artifacts 


The maximum normalized recovery error comparison 
between the proposed method and the ALO method [3] for 
different compression ratios for ECG signals is shown in Table 
I. Identical settings were maintained for the algorithms for both 
methods. The simulations were carried out in MATLAB 
R2014b environment with Intel i5-8250U CPU @1.6GHz and 
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8 GB RAM. ECG signal with motion artifacts and recovered 
signal using the proposed scheme for a compression ratio of 
50% is shown in Fig.8. 


From the results, it is clearly established that with the 
proposed CS scheme, ECG signals can be compressively 
acquired (with M much lower than N) and still be efficiently 
recovered from these compressive measurements without much 
losses. More significantly, from the results of ECG recovery for 
arrhythmia and motion artifacts cases where the ECG patterns 
are irregular or with significant pattern variations, the 
robustness of the proposed scheme is established. 


V. CONCLUSIONS 


A compressed sensing and recovery method is proposed for 
ECG signal using simple binary measurement matrix and 
inverted probability density function of Laplace Distribution as 
constraint function used along with traditional Gradient decent 
method. Normal ECG samples as well as anomalous cases like 
arrhythmia and ECG signals with motion artifacts from 
Physiobank ATM and MIT BIH databases were used to 
extensively evaluate the proposed scheme. The results were 
compared with the ALO method [3] proposed in literature for 
different compression ratios. The proposed work is found to 
yield improved results in compressively sensing and sparse 
recovery of ECG signals with good recovery characteristics. 
With lower computational resource requirements like lower 
storage requirement for the measurement matrix and no 
specific optimization requirement for measurement matrix, the 
proposed method is very suitable for real time CS applications. 
The proposed method can be advantageously used for efficient 
acquisition and sparse recovery of ECG signals in bandwidth 
constrained applications like manned spaceflights. 
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